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ABSTRACT 


In recent years, spacecraft have become increasingly flexible. The design 
requirements for the Bifocal Relay Mirror spacecraft include controlling jitter at 
the nanoradian level. Typically, tight pointing requirements require high 
structural stiffness, at the cost of increasing the on-orbit mass. To accomplish 
this, while minimizing the mass of the spacecraft, the structure will have some 
inherent flexibility. These flexible modes will interact with the pointing control, 
hence affecting the payload performance. The compensator design conducted in 
this thesis achieves order of magnitude improvements in controlling the rate 
error, hence jitter. This thesis starts with a rigid body dynamic model, and 
develops a flexible body dynamic model. Once the model is developed, the 
structure-controls interaction is discussed. Finally, compensators are applied to 
the rigid body controller to mitigate the performance losses present in the flexible 
body system. Through classical second-order compensators, the angular rate 
error was decreased by a factor of ten. Nonminimum phase notch filters and 
phase lag filters were used. Ultimately, the phase lag filters provided the best 


performance. 
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I. INTRODUCTION 


Recently, spacecraft designs have become more flexible in comparison to 
the legacy satellites, which are more rigid. There are also designs that involve 
multibody motion, such as a moving payload on a slewing central bus. The 
Bifocal Relay Mirror combines these situations for the purpose of relaying a laser 
from a ground or air source to space relay, to a ground or air target or receiver. 
The target may be cooperative, in the case of optical communications or other 
non-force application purposes, or uncooperative in the case of high-energy 


applications. 


Increasing the precision of the acquisition, tracking and pointing sequence 
demands that multibody dynamics and flexible dynamic effects be combined and 
addressed from a control design perspective. This thesis presents the rigid body 
and _ flexible body dynamics, and compares the uncompensated and 
compensated performance of the spacecraft using a flexible plant. Classical 


compensator designs are reviewed and applied. 


The design requirements for the Bifocal Relay Mirror spacecraft include 
controlling jitter at the nanoradian level. Typically, tight pointing requirements 
require high structural stiffness, at the cost of increasing the on-orbit mass. To 
accomplish this, while minimizing the mass of the spacecraft, the structure will 
have some inherent flexibility. These flexible modes will interact with the pointing 
control, hence affecting the payload performance. The compensator design 
conducted in this thesis achieves order of magnitude improvements in controlling 


the rate error, hence jitter. 


Through classical second-order compensators, the angular rate error was 
decreased by a factor of ten. Nonminimum phase notch filters and phase lag 


filters were used. Ultimately, the phase lag filters yielded the best performance. 
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ll. BIFOCAL RELAY MIRROR RIGID BODY EQUATIONS OF 
MOTION 


This section presents the derivations from Appendix B of the Bifocal Relay 
Mirror Technology Development paper submitted by the Naval Postgraduate 
School in February 2002 for the National Reconnaissance Office Directors 
Innovation Initiative competition. It provides the derivations of the equations of 
motion for the bifocal relay mirror spacecraft. The spacecraft sketch and relative 
size are shown in Figure 1. It consists of two bodies, transmit telescope and 
receive telescope. The receive telescope rotates with respect to the transmit 
telescope about a single axis. Figure 2 depicts a typical concept of operational 
employment for the Bifocal Relay Mirror Spacecraft. The receive telescope 
acquires and tracks the source laser, while the transmit telescope tracks the 
target. The motion of the fast steering mirrors accomplished internal alignment 
and fine pointing of the beam on the target. 
A. COORDINATE SYSTEM DEFINITION 

The coordinate frame descriptions are with respect to Figure 3. The 
center of mass of the receive telescope is on the rotation axis. Therefore, the 
center of mass of the system is fixed during the relative motion of the receive 


telescope. 










































































Figure 1 Bifocal Relay Mirror Spacecraft sketch (left) and relative size (right) 
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Figure 2. Concept of Operational Employment for Bifocal Relay Mirror Spacecraft 
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Figure 3 Bifocal Relay Mirror Spacecraft 
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The rigid body equations of motion for the transmitter and receiver are 
with respect to the inertial frame. Additionally, the dynamic coupling effects from 
the fast steering mirrors are neglected due to the small mass relative to the entire 
structure. 

The coordinate system 7,,7,,7, is fixed in the transmit telescope with 7, 
axis parallel to rotation axis of the receive telescope, 7, as telescope axis and 7, 
is normal to z, and rt, such that the 7,,7,,7, coordinate system is right-handed 
and mutually orthogonal frame. The origin of the coordinate system is at center 
of mass of the transmit telescope. 

The coordinate system /,, 0,, 0; is fixed in the receive system with 9, 
axis as the rotation axis, parallel to 7, axis, 6, as telescope and £/, is normal of 
p, and p, such that £,, ,,,' coordinate system is right-handed mutually 


orthogonal frame with origin at the center of mass of the receive telescope. The 


coordinate system /,, 0,, 03; is obtained from the coordinate system 7,,7,,7, by 
rotation a about x-axis. 

The relative motion of the transmit fast steering mirror with respect to 
transmit telescope is represented by rotation about 7, axis and £) rotation 7, 


axis. The relative motion of the receive fast steering mirror with respect to 


receive telescope is represented by £? rotation about 6, axis and #? rotation 
about , axis. The fast steering mirror rotations are assumed to be small. The 


equations of motion of the system are written in the coordinate frame ,, 7, 7; . 


B. EQUATIONS OF MOTION 
In a general case, rotational equations of motion for a given body about an 


arbitrary point P are given by 


Equation 1: Moment and Momentum Balance for Rigid Body 


M,=H,-p, x mf, 


p p 
where 


M. = total sum of external forces about P 


H,= angular momentum of the body about P 
p, = total vector from P to center of mass of the system 


m= _ total mass of the body 


r,= velocity of center of mass of the body 


If the point P is the center of mass, thenp, =0, thenM=H|,, with N 
being used to indicate the inertial reference frame. The angular momentum of 
the system, Hs can be written as follows: 

Equation 2: System Momentum Balance 


eee ee 


rel 


where 


H, = total angular momentum of the system. 


H-= total angular momentum by neglecting the contribution by the 


relative motion of the receive telescope and reaction wheels with respect to 


transmit telescope 


H_ = angular momentum due to relative motion of the receive telescope 


rel 


H., = angular momentum due to relative motion of the reaction wheels. 


To continue the momentum derivation, we must determine the inertia 


properties of the system. 


a. Inertia 
Let /, be the inertia matrix of the transmit telescope about its center of 


mass in coordinate frame z,, 7, 7, . 


i ae le 
belah sf oak, 
aly =f is 


Let the vector from center of mass of the transmit telescope to the center 


of mass of the spacecraft is given by 
pet = X.0 tyr) +z, 4 


Let I’ be the inertia matrix of the receive telescope about its center of 


mass in coordinate frame /,, 0,, P;. It is given by 


Ad Te 

x'x x'y x'Z 

ye | Be My le 
le TE TE, 

z'x zy z'z 


Next we transform L to z,7,,7,frame. The transformation matrix from 
Ts ia iz to Pp Pas Ps is given by 
1 0 0 
°C’ =|0 cos(a)_ sin(a) 
0 -sin(@) cos(a) 
The transformation matrix ‘C’ from /,, 6,, 6; tO 7,7,,7, is transpose 
of °C*. I’in z,,7,,7, system, /,, is given by I, = °C?’ °C" 
Let the vector from center of mass of the receive telescope to the center of 
mass of the system is given by 

















re 





| 
=X ,P,+Y,Po+Z, 3 


The inertia matrix of the spacecraft about its center of mass in coordinate 


frame 7,,7,,7, is given by 


2 2 2 20 = 
y % Zz. —X.Y, —X,Z, Yp = Zp XY p Xf p 
7 2 2 2 2 
lp=it+ is PN SCY) SAG V2 >| RIN | =X YO eX HV 2: 
= es 2 2 2 2 
XZ, YZ, x, - ve —X 2, —YV p2p x, + v3 
2. Angular Velocities 


The inertial angular velocity of the transmit telescope 
iso =|0, @, @,| ". The relative angular velocity of the receive telescope 
with respect to transmit telescope is given by’é°=[a@ 0 0] "The inertial 


. . pn =>T 5 eo . T 
angular velocity of the receive telescope “6” = "6" +'G’ =|@,+a a, o,| . 


3. Angular Momentum 
The terms of Equation 2 are further defined below. The angular 


momentum without receive telescope relative motion is given byH = 1,“@°. The 


relative angular momentum is given byH,,, = l'o°. 


Through substitution, we have 

Equation 3: Refined System Momentum Balance 

H, =1,%@" +1," +Hy 

4. Kinematical Relationships 

Using Equation 1, the equation of motion of the spacecraft is given by 


Equation 4: Rigid Body Moment & Momentum Defined 


Nor ry 
—H,+°@° xH, 


—H, =rate of change in inertial frame 


T 


oi, = rate of change in transmit telescope frame 


It should be noted that since /, is function of a, it is therefore time 


dependent. 
a. Transmitter Kinematics 
The momentum balance of the transmitter is as _ given 


as: H’ = Lo high The total moment applied to the transmitter is given by: 
Tay “iO Nec: a 


M.=M,+M, = a H’ = ale & thw) 


The subscripts C and G represent Control and Gimbal, 
respectively. The gimbal torque is due to the relative motion of the two 
telescopes. The receiver momentum is given 


by“H? =1,%@? ="C?l,°C'(“@" +6") = °C’, °C*(“&' +a). The moment applied 


to the receiver is given by: 


Mig "ad 
vy v Noy T t(N oer = 
M, =-M, => dt H?* — ae C*T, LG ( O + @)) 
j d T € =¢ = ao es T oT = 
eee atk Cl °C (NG + &)) + "BH x CPI, °C (NG +) 


=(*C?1,°C* + C71, °C*)(N@" +) + °C*1, °C" (“O" + a) 
+ NOT x CPL PC (NG +a) 
System-level open-loop control input is achieved through system- 


level dynamic equilibrium, as follows. 


Equation 5: Open-Loop Control Torque 
M. =-M, +(M, +M,) =M, +M. 
= {(7CP, CF + CPE, PC") NB +.@) + CPI, °C" +) 
+ NB x CPL PC (NO + &)} + {INO + NO" x (INO + hy} 
al (NG +a) +1, NO" + NO x (lO + hy + °C? °C'R)+ 'C’L, °C 
In later chapters, it will become necessary to have the acceleration 


of the telescopes in order to determine the effects of flexibility. Given that 


attitude and angular rates are observable with most sensors, the above equation 


may be solved for’ . 


Equation 6: Angular Acceleration of Transmitter 


%O = Iron, (Mo — NO x (lrora, “@" + My + °C?L,°C'G)—1 (NG +&)— C71, °C') 


b. Receiver Kinematics 
The inertial angular acceleration of the receiver’ is 
simply “a? = “a +a. 
Equation 7: Angular Acceleration of Receiver 


‘or? =" +a 


QO 
= Trorae (Me — “Ox (Iporay SO" + Ay + °C?P1, °C’) —-1 (“+ &)— °C’, °C'R) +4 
5. Wheel Control Laws 


Reaction wheels are typically controlled by a proportional-derivative (PD) 
controller of the formHy, =2 k,QcQae+Kig@e- This control law is developed in 
Wie. Where q,-is the ith element of the error quaternion and @, is the error 


between commanded and measured angular rate about the ith axis. The PD 
controller used in this work is of the same form, but with the parameters defined 
differently. The feedback control law used in this work is shown in Equation 8. 


Equation 8: Proportional Derivative Controller 


Ay =K We Me + Te) 
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Where k, is the controller gain and z,is the controller time constant. This 
creates a dynamic situation of the form/,@ + k,t,0,, + k,9,-9,- =O. The bandwidth 
of the controller goes as the square root of k,. Hence, increasing the bandwidth 


by a factor of two requires increasing the gain by a factor of four. 

6. Commanded Trajectory 

The commanded trajectory for this control system is based upon initial 
conditions for the joint angle motion as well as the transmitter acceleration and 
rates. The initial quaternion for the transmit telescope is given by q,(t,). The 


initial transmit telescope inertial angular velocity is given by “a‘(t,). The initial 


joint angle is given by a(t,). The initial joint angle is given by a(t,). 


Equation 9: Initial Conditions for Simulations 


0 
q, (to) = : 
1 
0 
N@' (ty) =| 0 | rad/sec 
0 
a(t,)=0 rad 


a(t,)=0 rad/sec 


The commanded attitude profile is based upon a typical engagement 
scenario. The inertial acceleration for the transmit telescope and the relative 
acceleration for the receive telescope is taken from the existing simulation. The 


values are described in Equation 10. 


Equation 10: Commanded Acceleration Profile for Simulations 


0.308 sin(0.0079f) 
‘@'(t)=| 7.896sin(0.0126t) | x10° rad/sec? 
~1.097 sin(0.0105f) 


di(t) = 1.1484 sin(0.0209t)x10~ rad/sec? 
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lll. MULTI-BODY FLEXIBLE EQUATIONS OF MOTION 


A. FLEXIBLE EQUATIONS OF MOTION — GENERAL CASE 
Given: flexible body, with a cantilevered attachment at point A, the 


equations of motion may be stated as follows: 
Equation 11: Cantilevered Dynamics 
Ge a hott la eo tte| 
M,, Map U, Ky, Kaa U, F, 
Where 
| = Internal Nodes 
A = Attachment Point 


M&K = Mass & Stiffness Matrices from cantilevered condition in finite 


element analysis program 
Equation 12: Transfer Nodes to Modal Coordinates 
vetf UL afPe Pal) % 
U, O 8 6/ |[U, 


Where q, are modal coordinates at the attachment and ®, is normalized 
for unity modal masses, i.e.®.'M,®, =/. The terms modal coordinates and 


generalized coordinates are interchangeable for the purposes used in this 
document. 


Z 
Substitute Equation 12 into Equation 11 and premultiply by Ee iI 
R 


get: 


Equation 13: Dynamics Expressed in Modal and Attachment 
Coordinates 


wr weet lS ollettte| 


Where, 
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Mae = 0,'M,®p 0.'M, 
Map = D,'M Pp 5 M,®p = O,'M ia 5 Maa 


o* = 0.'K,®, 
ia 0,'F, 
F,=0,'F,+F, 


Equation 13 is general in nature, and specific to the coordinates of the 
flexible body. The next series of equation transforms Equation 13 from 
coordinates at the cantilevered appendage to a common coordinate system with 
a given location and orientation. The selected coordinate system is defined as 
the axes defined by the transmit telescope coordinate system, located at the 
system center of mass, B. This basis is called {B}, as it is the body coordinate 
system use. The transmitter basis and the body basis differ in origin, but not 
orientation. The radius vectors from the B to the attachment points of the 


transmitter and receiver are given by: 


Equation 14: Center of Mass Locations 


foe xr tyr, +z.¢ 




















poe =X P:+¥,P>tZ, 3 
The attachment coordinates are expanded to divide between transmitter 
and receiver coordinates, superscripted by t and p respectively. Translation and 


rotation are indicated by lower case t and r respectively. For each flexible 
T 

body{U,}={U, Uj} . Hence, the translational and rotational attachment 

coordinates may be expressed in terms of the body coordinates as follows. 


Equation 15: Coordinate Transformations 


Ui ae Ie) Ie )-nm 


where 
0 Za Va 
[R|=|-z, 90 x, 
Ya —X 0 
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! 0 
Substitute Equation 15 into Equation 13 and premultiply by F | to 


obtain 


Equation 16: EOMs in Modal and Body Coordinates 


ie MES I-C2 
Mics. Maro U 0 0 U Fro 


B. FLEXIBLE EQUATIONS OF MOTION — TRANSMIT TELESCOPE 
Flexible Body Equations of Motion for Transmit Telescope (Note: A 


replaced by t) 


Coordinate Transfer: 
{T} = "C°{B} 

iz b, 

7, =4b, 

T3 b, 


‘OP =] 
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Equation 17: Transformation and Mass Matrices, Transmit Telescope 














400 0° 2. Hy, 
010 -z, O x, 
T]- 0.6 Aye xr 6 |! ‘ 
EO OO) ir Ge: HOF |e 4 
0000 1 0 
10000 0 1] 
—[(m 0 0 0 0 0 
0m 0 00 0 
O. 0. ae 20" 0 20.) ie, <0 
IMel=|q 9 9 ae oe |" | 
Or 0) 20. Jeo Gece 
iO, Oe 10! A Ie) 
[Meo] = 77M r.-| Me, Le | 
tRO T tRo et R! : R'M_R 4/7 


C. FLEXIBLE EQUATIONS OF MOTION — RECEIVE TELESCOPE 
Flexible Body Equations of Motion for Receive Telescope (Note: A 


replaced by p) 


Coordinate Transfer: 
{R} = °C*{B} 
ye} 1 0 0 b, 
Pop =|90 cos(a)  sin(@) |,b, 
Ps 0 -sin(a@) cos(a)||b, 
1 0 0 
’C® =|0 cos(a)_ sin(a) 
0 -sin(a@) cos(a) 
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Equation 18: Transformation and Mass Matrices, Receive Telescope 


EGY ACER 
ined 








1 0 0 0 Z, =V, 
0 cos(a) sin(a) y,sin(a)—Z,cos(a) -x,sin(a) x,cos(a) 
2 0 -sin(a) cos(a) y,cos(a)+z,sin(a) -x,cos(a) —x,sin(@) 
0 0 0 1 0 0 
0 0 0 0 cos(a) sin(a) 
| 0 0 0 0 —sin(a@) cos(a) 
[m, 0 0 0 0 O| 
Oo m, 0 0 0 O 
[ M 7 0. O.am- 0 0-0 7 M, 0 
7 OO) 402 fer IE. Ae IGS aie 
Or OF Oe Te Ie Te 
Oy 10" Oo. ie Ae IE 








[Mino |=T)M,2l, = 


p pR’ p 





eCPM PC® °C*TM °C8R 
p p p 
REPCFTM,°C® RT PCFM°C8R, + °CH I a 


Now, we have Maro for transmit and receive, we need Maco for transmit 
and receive. See Equation 16. These rigid-elastic coupling matrices are 
available from most of the widely used finite element analysis programs, such as 
NASTRAN or I-DEAS. 
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IV. STRUCTURE-CONTROLS INTERACTION 


In order to control the motion of flexible structures consisting of rigid body 
and elastic motion, the coupling between rigid body motion and elastic motion is 
important. Controlling the rigid body attitude is one of the primary objectives. 
Initially, the objective of achieving fine pointing may drive the control design to a 
simple, high gain controller. That is, retain the commonly used PD control 
architecture with higher proportional gain. This provides more control authority 
based to control attitude errors. If the designer increases the proportional gain 
by a factor of four, the control bandwidth doubles. This was discussed when 
Equation 8 was presented. If this high gain controller is within the structural 
frequency, the structural modes will get excited. Given the rigid-body model, the 
gain would be set in order to achieve the design-specific pointing accuracy. 
Ultimately, that is how the PD gains would be adjusted to achieve high pointing 
accuracy. The question is: how do we modify the controller to behave favorably 
around the structural frequencies? 

A. COMPENSATOR DESIGN 

Compensators are applied after the rigid body control is applied. The 
purpose of the compensator is to add stability to the system. The stability is with 
respect to gain, phase, or a combination. Classical gain-phase stabilization 
techniques are of primary use to single-input single-output (SISO) systems. A 
rigid or flexible body is said to be gain-phase stabilized if it is closed loop stable 
with finite gain and phase margins. These classical techniques will be applied to 
this multi-input multi-output (MIMO) system [Ref Wie]. 

1. Gain Stabilization 

Gain stabilization of a flexible body attenuates the gain at desired 
frequencies, to provide closed loop stability regardless of phase uncertainty. A 
gain stabilized mode has a finite gain margin, but is closed loop stable regardless 


of the phase uncertainty. 
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2. Phase Stabilization 

Phase stabilization of a flexible body provides phase characteristics at the 
desired frequency to obtain closed loop damping higher than the passive 
damping of the mode. A phase stabilized mode has a finite phase margin, but is 
closed loop stable regardless of the loop gain uncertainty. Since the dynamic 
model has no passive damping, phase stabilization is required for this 
application. 
B. COMMON CLASSICAL COMPENSATOR DESIGN 

Wie presents a comprehensive discussion compensator design in his 
book. One fundamental concept is that gain stabilization may only be achieved if 
the corresponding mode has some passive damping. Since no modal damping 
is modeled, phase stabilization techniques must be employed. The following 
sections highlight some common compensator designs. 


Control theory has made significant advances in the past several decades. 
Application of these theories has been difficult, from a practicality perspective. 
More often than not, classical techniques are applied. Performance is 
“optimized” through trial and error synthesis using frequency domain methods or 
actual performance in situ. Due to the complexity of the dynamic model, 
performance will be compared based on the performance of the plant under 
various compensation techniques. For lower performance requirements, simple 
first order filters are sufficient. The effect of the first order filter is to change the 
controller bandwidth. Typically, low pass filters are used for this type of an 
application, since they will eliminate the high frequency content of the controller, 
hence not exciting the structural modes. However, for higher performance 
requirements, such as placing the control bandwidth close to or within the 
structural modes, higher order filters become necessary. 


Wie presents the generalized second order filter, using the s-domain. He 
presents it as a natural extension of the classical notch or phase lead/lag filter, 
and is based on pole-zero patterns. The following sections employ Wie’s 
descriptions for various compensation methods. The filter is of the form: 
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Equation 19: Generalized Second-Order Filter, Ref Wie 





ec ee 
Oo (00) 
H(s)=— z 
26.8 
2 Sp +1 
Op , 
1. Minimum-Phase Lead Filter 


A phase-lead filter utilizes Equation 19, and assigns¢, =¢, =¢, >0. The 


maximum phase lead, obtained atw, = ,/a,@, , is determined by: 


Equation 20: Maximum Phase from a Phase Lead Filter, Ref Wie 


(2¢, Jo, /o, y -(o,/0, -1) 
(2¢. Jo, /o, i +(o,/0, ay 


The gain increases or decreases at high frequencies, and 


, at O, = .J0,0, 





Dice = cos | 


becomes K,, = 40log,,(@,/@,) dB. For small¢,, when filter poles and zeros are 
close to the imaginary axis, the effective lag region lies between, and a, , and 
the maximum phase shift approaches +180°. In practice, ¢, is often selected 
greater thanc,. This places the filter poles sufficiently far to the left of the 
imaginary axis. The use of a phase-lag filter with a large w,/a@, ratio greater 


than two should be avoided from a practical viewpoint. The phase lag filter is 


shown in Figure 4. 
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Bode Diagram 
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Figure 4 Bode Plot of Minimum-Phase Lag Filter 
2. Minimum-Phase Notch Filter 


Fora, =@,, a notch filter is obtained. The minimum gain of the filter is 
obtained at w, =@,=a@, asK,, =20log,,(¢,/¢,) dB. Both phase lead and lag 
occur neara@,. For the notch filter, the maximum phase lag and lead occur at a, 


and @,, respectively, where: 


Equation 21: Location of Maximum Phase Lead or Lag for Notch 
Filter, Ref Wie 


i = 20.6, +1-,|(2¢,6, +1) -1 
(4a) 


Cc 


2 = [2¢.6,41+, (26,6, +1) —1 
O 


Cc 
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Because 7 and % depend only on the product ¢,¢,,, the filter damping 
QO, QO, 


ratios determine the effective notch region. Typical notch filter damping ratios 


arec, =0,and ¢, =1. The minimum phase notch filter is shown in Figure 5. 
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Figure 5 Bode Plot of Minimum Phase Notch Filter 


3. Nonminimum-Phase Notch Filter 


Fora, =,, andg, >|¢,|(¢, <0), a nonminimum-phase notch filter is 


realized. The applicability of this filter is when passive damping does not exist 
(the ideal case). The conventional notch filter cannot be used in this case. As 
the pole-zero combination is placed further from the imaginary axis, the 
robustness of the filter is enhanced, while the stability margin is reduced. The 


nonminimum phase notch filter is shown in Figure 6. 
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Bode Diagram 
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Figure 6 Bode Plot of Noninimum Phase Notch Filter 


C. CLASSICAL COMPENSATOR DESIGN 

Wie presents an algorithm for designing compensators with generalized 
filters. The basic idea is to synthesize the compensator for one mode at a time. 
The design algorithm is summarized below: 

1: Select Control Bandwidth 

The driving control bandwidth is one of the key parameters of the design. 
This may be driven by performance, noise sensitivity, control authority, etc. 

2. Rigid Body Mode Compensation 

Rigid body control is typically achieved through a proportional-derivative 
(PD) feedback controller. The format of the control law is given 
by:u =—K(1+Ts)y =-(K,+K,s)y . 
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3. Flexible Mode Compensation 

The types and degrees of the closed loop behavior depend on the relative 
spectral separation of the modes. Determine if simple gain stabilization of all the 
flexible modes are possible. If not, determine the needed lead or lag angles foe 
each mode. Synthesize the appropriate structural filter for each destabilizing 
mode one by one, using the various second-order filters. 

4. Design Iteration 

Repeat the process to compromise the interactions between each 
compensator. A few iterations using a software package, or analytical model, will 


result in a more streamlined design process. 
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V. DEVELOPMENT AND VALIDATION ANALYTICAL MODEL 


A. FLEXIBILITY SIMULATION ARCHITECTURE 

Now that the flexible telescopes are modeled dynamically, a method of 
simulation must be developed. The approach used is the feedback torque effect. 
This method is presented in simple terms in Figure 7. In essence, the rigid body 
dynamics are unchanged. However, the state vector determines the base 
excitation (acceleration) which creates a feedback torque. This feedback torque 


is combined with the control torque to affect the rigid body dynamics. 





Rigid Body 
Dynamics 


Flexible Body 
Dynamics 


Figure 7 Flexibility Model, After Ref Agrawal 





In Figure 7, T represents the command input (Torque), and 6) represents 
the state vector. Using Equation 16 as a reference, the general equation that we 
are solving is of the form Ax + Bx =F, where x = [q, U]'. The most practical and 
mathematically stable method of solving this system of equations is to determine 
the rigid body state and work backwards to determine the reaction torque which 
results from the flexibility. This is shown in the next series of equations. If we 


start with Equation 16 as a base equation, the known values are the mass and 
stiffness matrices, as well as U and q,- We may restate Equation 16 as two 


differential algebraic equations, as in Equation 22. 
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Equation 22: Application of Equation 16 for Determining Simulation 
Algorithm 


IG, +MygoU + aq, =0 
Maco Ga + MagoU = T, 


Through reorganization, Equation 22 may be restated as 


Equation 23: Determining Reaction Torque 


Va = —M geoU rs OF, 


reaction — —(Maco Ga a MypoU) 

At first, it isn’t intuitive why the second equation has the negative sign. 
Without that negative sign, the system is unstable. With the negative sign, the 
system behaves as it should. In plain verbiage, the negative sign is there 
because it is a reaction torque, the behavior of which opposes that of an applied 
torque. The generalized coordinates are initialized at zero when the simulation 
begins. The stable algorithm for determining the reaction torque is shown in 


depicted in Figure 8. 








R-E Coupling 
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1 
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Figure 8 Determining Feedback Torque from Rigid Body Motion 


B. DYNAMICAL SIMULATION DEVELOPMENT 

The base acceleration is determined using Equation 6 and Equation 7 for 
the transmitter and receiver, respectively. The block diagram in Figure 9 shows 
how to determine the inertial acceleration of each body from the observed states, 
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using Equation 6 and Equation 7. The block diagram in Figure 10 determines the 


feedback torque due to flexibility, using Equation 23. 
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Figure 9 SIMULINK Model for Determining Rigid Body Acceleration of Transmitter 
and Receiver, using Equation 22 and Equation 23 
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Figure 10 SIMULINK Model for Determining Torque Due to Acceleration of Flexible 
Bodies, using Equation 23 
C. VALIDATION OF DYNAMICAL MODEL 
1. Developing Flexible Bodies 
The MSC software suites of PATRAN and NASTRAN were used for 
modeling the geometry and determining the modes of each telescope. The 
objective was to have a first mode for the transmitter near 10 Hz, with the 


receiver being about slightly higher, due to less mass and the same stiffness. 
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This requirement was changed to around 1 Hz for the purposes of propagating 
an answer. Lower modal frequencies take less time to propagate. 
a. PATRAN Telescope Geometry 


| Transmitter | Receiver_| 
Inner Diameter (m) 
Outer Diameter (m) 
Length (m) 


Mass (kg) 
Density (kg/m’3) 
Young's Modulus (Pa) 


Figure 11 Flexible Body Parameters from Geometry in PATRAN 





b. NASTRAN Normal Modes 

The resultant modes from NASTRAN are shown in Figure 12. 
During simulation, the impacts of flexibility were impractical for such high 
frequencies. Therefore, the modes were modified with a flexibility factor. It was 
determined that reducing the modes by a factor of ten resulted in a manageable 


simulation. The modes used for simulation are shown in Figure 13. 


| 6 | 39.58 | 60.35 _| 
| 8 | 49.67 | 75.75 
| 9 | 54.92 | 79.18 | 





Figure 12 Receiver and Transmitter Modes from NASTRAN (Hz) 
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| 6 | 396 | 6.04 | 
| 8 | 407 | 7.58 
| 9 | 519 | 7.92 | 





Figure 13 Receiver and Transmitter Modes Used in Simulation (Hz) 


Cc. Modeling Rigid-Elastic Coupling 

The rigid-elastic coupling matrix was not taken from a finite element 
analysis program. Instead was chosen arbitrarily, from visual inspection of the 
displayed modes in the simulation mode of PATRAN. However, the method 
used could easily be substituted by NASTRAN output. The arbitrary method 
used involved selecting a unit vector for each mode considered, and then 
multiplying this by a coupling factor. The unit vectors represent the amount of 
the coupling that is occurring between an axis of rotation and a generalized 
coordinate. The first two modes are bending modes about the x and y-axis 
respectively. The third mode is a twisting mode about the z-axis. The fourth and 
fifth modes are combinations of these effects. The rigid-elastic coupling matrix 
was adjusted such that the maximum coupling occurred from the first mode, and 
the minimum coupling occurred from the last mode considered. The coupling 
factors varied from a maximum of 10% to a minimum of 1% of the maximum 


principal moment of inertia, and are shown in Figure 14. 


Coupling Factor 
7000 
5425 


3850 
2275 





Figure 14 Rigid-Elastic Coupling Factors for Simulation 
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The rigid-elastic coupling matrix is presented in the following, where 
Mago = Moco = M.c0- 


Myrco = CouplingFactorMatrix * Coupling Vectors 




















[7000 0 0 0 0 1 0 0 
O 5425 0O 0 0 0 1 0 
=| O 0 3850 0 0 0 0 1 
0 0 0 2275 O || 0.7071 0 0.7071 
| O 0 0 O 700}|0.7071 0.7071 0 
| 7000 0 0 | 
0 5425 0 
= 0 0 3850 |kg—m? 
1608.67 0 1608.67 
[494.975 494.975 0 
2. Observing System Free-Free Behavior 


Nominally, all generalized coordinates are assigned an initial condition of 


zero. In order to demonstrate the uncontrolled dynamics, generalized 
coordinates are given random initial conditions, and no actuation is applied. The 
telescope joint angle motion and spacecraft attitude actuators (reaction wheels) 
are set to zero input. This is equivalent to stretching and compressing various 
springs in an undamped multi-body spring-mass system and letting them go. 


The purpose is to observe the free-free behavior of the system. 


Figure 15 was created with random initial conditions for the generalized 
coordinates. Additionally, the initial angular rate of the transmitter was [0,0,0] 
rad/sec. There was no commanded angular acceleration of the transmitter or 
receiver, and there was no initial angular momentum. The purpose of this figure 
is to demonstrate the complex frequency content of the response. The order of 
magnitude of the difference in angular rates is not to be compared to future 
graphs. Since the initial conditions are random, any such comparison would be 


meaningless. 
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x10° Difference in Angular Rates 


T T T T 














Figure 15 Uncommanded Free-Free Behavior, Angular Rates 


3. Eigenvalues of the Free-Free System 

Fundamentally, the free-free behavior of the system is dependent on the 
joint angle. This assumes that the rate of motion has no effects on the free-free 
modal behavior. Since the input has no effect on the modes, the base equation 


is shown in. 


Equation 24: Free-Free Dynamics, No External Forces 


/ 0 Maco || 9: o, 0 O}}4q, 0 
0 | M seo q, 0 o, 0 qV,¢= 0 
M ee M “eas Maro || U 0 0 O;;|U 0 


When this is put into the eigenvalue format, we have Equation 25. 
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Equation 25: Free-Free Dynamics, Eigenvalue Solution 
det {271-M'kK} =0 


where 
0 M seo 
M=| 0O / M seo 
Mico’ Mico Maro 
a> > @ 
K=]| 0 o,” 0 
0 0 O 


A simulation was created to determine the eigenvalues of the system at a 
family of joint anglesa<[0 z]. This simulation solves Equation 25, and is 


shown in Figure 16. 
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Figure 16 SIMULINK Model for Determining Eigenvalues of the Free-Free System, 
using Equation 25 


The eigenvalues are sensitive to the coupling factors for the rigid-elastic 
coupling. As the coupling factors are reduced, the system free-free frequencies 
approach the modes of each single body. These free-free modes are what must 
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be compensated for in the control system. For the coupling factors used, the 
free-free system frequencies are shown in Figure 17. The column identified as 
“without coupling” uses coupling factors of zero. In all columns, the rigid body 
modes are omitted, as the frequencies are zero for these modes. 


Frequencies (rad/sec 


) 
List of Single 
Body Modes 


Figure 17 Free-Free Modes with and without Rigid-Elastic Coupling 





D. SENSITIVITY OF PLANT BEHAVIOR TO UNCOMPENSATED PD 
CONTROLLER BANDWIDTH 


This section shows the interaction between the PD controller bandwidth 
and the flexible plant, without any compensator applied. This is accomplished by 
turning off all compensators in the model, holding the time constant unchanged, 
and steadily increasing the controller gain. 


Fundamentally, designing a control law for a system with flexible bodies 
seeks high bandwidth control to achieve fine pointing, while not exciting flexible 
modes. When the original rigid body gains for the PD controller were applied to 
the flexible system, they were too low. The attitude and rate errors which 
occurred were due to a low controller gain, and not vibration. To achieve higher 
pointing accuracy, the controller gains were increased while the time constants 


were unchanged. 
The original rigid body controller had gains around 5000 and time constant 
around 0.3 sec. The time constant was adjusted to t = 0.01sec in order to lower 


the damping of the rigid body controller. Then, the gains were adjusted to find 
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the points of instability for various modes. The highest stable gain was 
determined, as well as the gain which caused instability for the first three modes. 
The adjusted values of K increase the bandwidth of the controller, and 


demonstrate structure-controls interaction. 


The figures that follow show the results of the uncompensated flexible 
body for increasing gains. There are three plots in each figure. The top plot is 
the difference in quaternions. This has no physical meaning, but does allow the 
user to see the rate at which error is occurring. The middle plot is the difference 
in angular rates: commanded rate minus actual rate. The bottom plot is the 
pointing error, in radians: angle subtended between the actual attitude and the 
commanded attitude. 


In Figure 21, the angular rate errors are on the order of 1e-6 rad/sec, and 
the pointing error continues to grow. This is an effect of a low gain. In Figure 22, 
Figure 23 and Figure 24, the angular rate errors continue to decrease, and the 
onset of angular error is delayed. Finally, in Figure 25 there is no angular error 
throughout the simulation. Now the gain is increased until the plant cannot be 
controlled. The rate error continues to decrease in Figure 26 and Figure 27. 
However, in Figure 28, the system instability is obvious. 


Through trial and error, the gains at which each axis want unstable were 
determined. Figure 29 shows the behavior at the highest stable gain, 50.3e6. At 
this gain, the angular rate error is on the order of 5e-10 rad/sec. Uncompensated 
gains above this level result in rate errors on the order of 5e-4 rad/sec, one 
million times higher than the highest stable gain. The first instability occurs at 
50.4e6, as shown in Figure 30. The second axis is unstable at 85.7e6, as shown 
in Figure 31. The third axis is unstable at 93.6e6, as shown in Figure 32. For all 
of these simulations, the initial pointing error was 10 nanoradians. That causes 
the initial oscillation. 
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Max Rate Error | Angle Error at 
at 20 seconds 20 seconds 
Gain (rad/sec) (rad 


5,000.00 5.E-06 3.E-05 


20,000.00 1.E-06 8.E-06 
80,000.00 2.E-06 


1,280,000.00 
5,120,000.00 
20,480,000.00 
[50,300,000.00] 5.E-10 | 0 


Figure 18 PD Controller Gain and Corresponding Errors for Flexible Spacecraft, 
Without Compensator 


Max Rate Error at] Angle Error at 
Compensator 20 seconds 20 seconds 
Gain Type Omega_z|Omega_p| Zeta z Zeta_p (rad/sec) (rad) 
2 


51,000,000|NMP Notch _ [first three modes -0.00001 5.00E-0 
60,000,000|NMP Notch _ |first three modes -0.00001 4.00E-07} ———*OJ 
55,000,000|NMP Notch _ [first three modes -0.00001 3.00E-07/| OO 


Figure 19 PD Controller Gain and Corresponding Errors for Flexible Spacecraft, 
Nonminimum Phase Notch Filter 


Max Rate Error at} Angle Error at 
Compensator 20 seconds 20 seconds 
Gain Type Omega_z}Omega_p| Zeta _z Zeta_p (rad/sec) (rad) 
ee ee a) 


| 60,000,000|PhaseLag | 14] 5] 4 O0E-10] 
| 86,000,000|PhaseLag | 14] 5] 38.0010] 
| 95,000,000|PhaseLag_ | 14] tT 5] 3. OOE-10] — 
|105,000,000|PhaseLag | 14] 5] 3. 00E-10] 
[115,000,000|PhaseLag | 14] 5] 200-10] 
1125,000,000|PhaseLag | 14] 5] 200-10] 
[140,000,000|PhaseLag | 14] 5] B0E-10] 
|160,000,000|PhaseLag_ | 14] SY] 5OE-10] 
|258,000,000|PhaseLag | 14] 5] 1 00E-10] 
|258,000,000|PhaseLag_ | 20] G3] 0.33] 1. 20E-10] 
|258,000,000|PhaseLag_ | 20] SG] 3] 0385] 20E-10] 
|258,000,000|PhaseLag_ | 20] G21 20E-10] 
|500,000,000/PhaseLag | 20] tS OOE-T1] 


Figure 20 PD Controller Gain and Corresponding Errors for Flexible Spacecraft, 
Phase Lag Filter 


320,000.00 3.E-07 
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Figure 21 Uncompensated Behavior, K = 5000 
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Figure 22 Uncompensated Behavior, K = 20,000 
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Figure 23 Uncompensated Behavior, K = 80,000 
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Figure 24 Uncompensated Behavior, K = 320,000 


40 


x 10° Difference in Quaternions 
































10 t ; t t 
BL 4 
Opesesee 
5 L 1 1 1 L L i 1 1 
0 2 4 6 8 10 12 14 16 18 20 
x 10° Difference in Angular Rates 
2F T T T T T 7 
tHi| | | 4 
Aft | 
olf VA ; 
| } 
Vey 
14 a 
L 1 1 1 L L 1 1 L i 
0 2 4 6 8 10 12 14 16 18 20 
Angle Error (rad) 
1 t ; t t 
0.5; 4 
0 
-0.5}- 4 
| L | | | L L | | | 
0 2 4 6 8 10 12 14 16 18 20 


— xaxis 
— y-axis 
— z-axis 


Figure 25 Uncompensated Behavior, K = 1,280,000 
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Figure 26 Uncompensated Behavior, K = 5,120,000 
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Figure 27 Uncompensated Behavior, K = 20,480,000 
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Figure 28 Uncompensated Behavior, K = 81,920,000 
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Figure 29 


Figure 30 
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Uncompensated Behavior, K = 50,300,000, the last stable gain 
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Uncompensated Behavior, K = 50,400,000, two axes unstable 
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Figure 31 Uncompensated Behavior, K = 85,700,000, onset of second axis unstable 
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Figure 32 Uncompensated Behavior, K = 93,600,000, onset of third axis unstable 
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VI. COMPENSATOR APPLICATION TO ANALYTICAL MODEL 


Based on the compensators described in Chapter IV, two applicable 
designs were selected for application to this control system. The Nonminimum- 
Phase Notch Filter and the Minimum-Phase Lag Filter were chosen for this 
undamped system. The conventional Minimum Phase Notch Filter is not 
applicable since there is no passive system damping for the flexible modes. The 
structural filters used in industry are of the Minimum-Phase Lead variant. This 
type of filter was attempted. However, no performance increases were realized. 
The lowest error achieved without a compensator occurred at K = 50.3e6. At this 
gain, the angle error is numerically “zero”. Although, it is not zero, it is below the 
threshold of which MATLAB is able to take the arccosine. The rate error at the 
end of the simulation was approximately 5e-10 rad/sec. The objective of the 
compensator is to reduce this rate error, thereby reducing the impact of structural 
vibrations on the pointing error budget. 

A. NONMINIMUM PHASE NOTCH FILTER 

The first method used was the nonminimum phase notch filter. The 
compensators were designed about each identified free-free frequency. Initially, 
one was applied. By the end, three compensators were applied, with the 
purpose of suppressing the first three modes. The best performance achieved 
was on the order of 3e-8 rad/sec. Unfortunately, this is worse than the 
uncompensated system before the onset of instability. The results are shown 


below. 
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Figure 33 Notch Filter, K = 51e6, ¢, = —0.00001, ¢, =1.0 
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Figure 34 Notch Filter, K = 60e6, ¢, = —0.00001, ¢, =0.7 
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Figure 35 Notch Filter, K = 55e6, ¢, = —-0.00001, ¢, =0.7 








B. MINIMUM PHASE LAG FILTER 

The second method used was the minimum phase notch filter. The 
compensators were designed about the identified free-free frequency. The best 
performance achieved was on the order of 5e-11 rad/sec. This is a factor of ten 


reduction in angular rate error. 


Initially, a highly damped filter was used. With ¢,=1 and ¢,=1.5, the gain 


was steadily increased from Figure 36 through Figure 45. The rate error in 
Figure 45 was reduced to 1e-10 rad/sec. At this point, with a gain of 258e6, the 


error has been reduced by a factor of five. 


After this, the filter was lightly damped. The zeros had damping ratios 
between 0.2 and 0.3. The poles had damping ratios 1.05 times that of the zeros. 


Different cases are presented in Figure 46 through Figure 49. The rate error in 
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Figure 49 was reduced to 5e-11 rad/sec. At this point, with a gain of 500e6, the 


error has been reduced by a factor of ten. 
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Figure 36 Phase Lag Filter, K=51e6, ¢€,=1, ¢,=1.5, w,=14, w,=6 
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Figure 37 Phase Lag Filter, K=60e6, ¢,=1, ¢,=1.5, 7,=14, w,=6 
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Figure 38 Phase Lag Filter, K=86e6, ¢,=1, ¢,=1.5, 7,=14, w,=6 
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Figure 39 Phase Lag Filter, K=95e6, ¢,=1, ¢,=1.5, ,=14, w,=6 
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Figure 40 Phase Lag Filter, K=105e6, ¢,=1, ¢,=1.5, @,=14, w,=6 
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Figure 41 Phase Lag Filter, K=115e6, ¢,=1, ¢,=1.5, @,=14, @,=6 
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Figure 42 Phase Lag Filter, K=125e6, ¢,=1, ¢,=1.5, @,=14, w,=6 
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Figure 43. Phase Lag Filter, K=140e6, ¢,=1, ¢,=1.5, 7,=14, w,=6 
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Figure 44 Phase Lag Filter, K=160e6, ¢,=1, ¢,=1.5, @,=14, w,=6 
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Figure 45 Phase Lag Filter, K=258e6, ¢,=1, ¢,=1.5, 7,=14, w,=6 
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Figure 46 Phase Lag Filter, K=258e6, ¢,=0.3, ¢,=0.33, w,=20, w,=6 
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Figure 47 Phase Lag Filter, K=258e6, ¢,=0.33, ¢,=0.35, w,=20, w,=6 
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Figure 48 Phase Lag Filter, K=258e6, ¢,=0.2, ¢,=0.21, w,=20, w,=6 
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Figure 49 Phase Lag Filter, K=500e6, ¢,=0.2, ¢,=0.21, w,=20, w,=6 
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Vil. CONCLUSION 


The design requirements for the Bifocal Relay Mirror spacecraft include 
controlling jitter at the nanoradian level. Typically, tight pointing requirements 
require high structural stiffness, at the cost of increasing the on-orbit mass. To 
accomplish this, while minimizing the mass of the spacecraft, the structure will 
have some inherent flexibility. These flexible modes will interact with the pointing 
control, hence affecting the payload performance. The compensator design 
conducted in this thesis achieved order of magnitude improvements in controlling 
the rate error, hence jitter. This thesis presented the rigid body and flexible body 
dynamics, and compares the uncompensated and compensated performance of 
the spacecraft using a flexible plant. Classical compensator designs were 
reviewed and applied. 


The plant in this work was nonlinear. No attempts were made to linearize 
the plant about a set of conditions, or for a given set of command inputs. The 
instabilities were found completely by trial and error. One alternate approach for 
determining the destabilizing gain would be to command a rotation about one of 
the principal axis. Even with a time variant parameter, this would provide a linear 
plant. This would allow the control designer to use classical controller design, 
such as root locus techniques, to determine the unstable gains. This thesis did 
no such linearization, as a primary task was to establish a baseline nonlinear 


simulation to which others could be compared. 


This thesis showed that classical, high gain PD controllers are appropriate 
for rigid body control. However, when the controlled plant includes flexible 
modes, uncompensated high gain PD control causes extreme deviation due to 
unmodeled dynamics. The compensators used in this formulation validate their 
application on flexible spacecraft with high pointing accuracy, such as the Bifocal 
Relay Mirror. The relative motion of the two telescopes is slow in this case. No 
effort was made to adjust the compensator based on the joint angle motion 
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characteristics. Classical compensator synthesis was a sufficiently accurate 
method for controlling the feedback torque due to flexibility. 


Through classical second-order compensators, the angular rate error was 
decreased by a factor of ten. Nonminimum phase notch filters and phase lag 


filters were used. Ultimately, the phase lag filters provided the best performance. 
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